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ABSTRACT 

Aims. LS 5039 is a High Mass X-ray Binary (HMXRB) comprising a compact object in an eccentric 3.9 day orbit around a 
massive 06.5V star. Observations at energies above 0.1 TeV (10 11 eV) by the High Energy Stereoscopic System (HESS) in 2004 
revealed that LS 5039 is a source of Very High Energy (VHE) 7-rays and hence, is able to accelerate particles to multi-TeV 
energies. Deeper observations by HESS were carried out in 2005 in an effort to probe further the high energy astrophysics 
taking place. In particular, we have searched for orbital modulation of the VHE 7-ray flux, which if detected, would yield new 
information about the complex variation in 7-ray absorption and production within X-ray binary systems. 
Methods. Observations at energies above 0.1 TeV (10 1 eV), were carried out with the High Energy Stereoscopic System (HESS) 
of Cherenkov Telescopes in 2005. A timing analysis was performed on the dataset employing the Lomb-Scargle and Normalised 
Rayleigh statistics, and orbital phase-resolved energy spectra were obtained. 

Results. The timing analysis reveals a highly significant (post-trial chance probability < 10 -15 ) peak in the TeV emission 
periodogram at a frequency matching that of the 3.9 day orbital motion of the compact object around the massive stellar com- 
panion. This is the first time in 7-ray astronomy that orbital modulation has been observed, and periodicity clearly established 
using ground-based 7-ray detectors. The 7-ray emission is largely confined to half of the orbit, peaking around the inferior 
conjunction epoch of the compact object. Around this epoch, there is also a hardening of the energy spectrum in the energy 
range between 0.2 TeV and a few TeV. 

Conclusions. The 7-ray flux vs. orbital phase profile suggests the presence of 7-ray absorption via pair production, which 
would imply that a large fraction of the 7-ray production region is situated within ~1 AU of the compact object. This source 
size constraint can be compared to the collimated outflows or jets observed in LS 5039 resolved down to scales of a few AU. 
The spectral hardening is however not explained exclusively by the absorption effect, indicating that other effects are present, 
perhaps related to the 7-ray production mechanism (s). If the 7-ray emission arises from accelerated electrons, the hardening 
may arise from variations with phase in the maximum electron energies, the dominant radiative mechanism, and/or the angular 
dependence in the inverse-Compton scattering cross-section. Overall, these results provide new insights into the competing 7-ray 
absorption and production processes in X-ray binaries. 



1. Introduction 

X-ray binaries (XRBs) comprise a compact object such 
as a neutron star or black hole in orbit with a compan- 
ion star. They are one of several types of astrophysical 
system that can provide a periodic environment for the 
acceleration of particles and subsequent production of ra- 
diation. Modulation of this radiation, linked to the orbital 
motion of the binary system, provides key insights into 
the nature and location of particle acceleration and emis- 
sion processes. While such modulation is often found in 
XRBs up to hard X-ray energies (Lewin I1995|, Wen et al. 
2006), until now it has not been established in any as- 
trophysical source at 7-ray energies. LS 5039 (distance 
d ~2.5 kpc) is a HMXRB comprising a compact ob- 
ject in a ~3.9 day orbit around a massive 06.5V star 
(Casares et al. I2005|) . Persistent radio outflows (observed 
with extension in the range 2 to ~1000 AU) are attributed 
to a mildly relativistic (v ~ 0.2c) jet (Paredes et al. 
2000, 2002 ), which would place LS 5039 in the microquasar 
class. Microquasars are Galactic, scaled-down versions of 
Active Galactic Nuclei (AGN) (Mirabel fT994]) . and are a 
sub-class of XRBs. The detection of radio (Marti et al. 
119981 Ribo et al. [1999) and X-ray (Bosch-Ramon 
et al. 2005) emission and their possible association 
with the MeV to GeV 7-ray sources GRO J1823-12 
(CollmarEOOal) and 3EG 1824-1514 (Paredes et al. 12000]) 
suggests the presence of multi-GeV particles. Observations 
in 2004 (Aharonian et al. I2005b|) (~11 hrs) with HESS 
established a new VHE 7-ray source, HESS J1826-148, 
within 30 arcsec of the radio position of LS 5039, re- 
vealing for the first time that LS 5039, and hence XRBs, 
are capable of multi-TeV (10 12 eV) particle acceleration. 
The limited statistics did not allow for detailed timing 
or variability analyses. We note that evidence for vari- 
ability at VHE 7-ray energies has recently been unveiled 
(Albert et al. [2006) in a similar type of binary system, 
LS I +61°303. Here we report on new, deeper HESS ob- 
servations of LS 5039 at TeV 7-ray energies, revealing that 
its VHE 7-ray emission is modulated by the orbital motion 
of the compact object around its massive stellar compan- 
ion. 

2. Observations 

The observations were taken with HESS (Aharonian 
et al. [2006b). an array of four identical Atmospheric 
Cherenkov Telescopes (ACT) located in the Southern 
Hemisphere (Namibia, 1800 m a.s.l), and is sensitive 
to 7-rays above .1 TeV . The 2004 HESS observations 
(Aharonian et al. I2005bp have been followed up with a 
deeper campaign in 2005. After data quality selection, the 
total dataset comprises 160 runs (or pointings) represent- 
ing 69.2 hours observations from both 2004 and 2005. Data 
were analysed, employing two separate calibration proce- 
dures (Aharonian et al. I2004|) and several background re- 

* Present address: School of Chemistry and Physics, 
University of Adelaide, 5005, Australia 



jection and direction reconstruction methods. The results 
presented here are based on the combination of a semi- 
analytical shower model and a parametrisation based on 
the moment method of Hillas to yield the combined like- 
lihood of the event being initiated by a 7-ray primary 
(de Naurois et al. l2003"l) . As we show later, a pure Hillas- 
based analysis, described in Aharonian et al. (2006b), also 
yielded consistent results. 

3. Results 

A total of 1960 7-ray events (with an excess significance 
above the background exceeding +40er) within 0.1° of the 
VLBA radio position of LS 5039 (Ribo et al. [2002)) were 
found. The best-fit position (in Galactic coordinates) is 
I = 16.879°, b = —1.285° with statistical and systematic 
uncertainties of ±12 and ±20 arcsec, respectively, which 
is consistent with the VLBA position within the la sta- 
tistical uncertainty. 

3.1. Timing Analysis 

A search for periodicity, by decomposing the runwise VHE 
7-ray flux at energies > 1 TeV into its frequency com- 
ponents, was carried out using the Lomb-Scargle Test 
fScargle [T98"2"|) . (Fig. [J) and Normalised Rayleigh Statistic 
(NRS) (de Jager I1994[) (Fig. [2]) which are appropriate 
for unevenly sampled datasets typical of those taken by 
HESS. The 2005 observations were taken over a wide range 
of zenith angles yielding a varying energy threshold in the 
range 0.2 to ~1 TeV. To reduce adverse affects of this 
varying threshold in our timing analysis, wc used all events 
and extracted the flux normalisation above 1 TeV assum- 
ing an average photon power-law index derived from all 
data (r = 2.23 for dN/dE ~ £T r ). Although as we see 
later the photon spectral index was found to vary within 
the orbital period, the average index assumption in this 
method contributes only a small error on the derived flux 
above 1 TeV. 

An obvious peak in the Lomb-Scargle periodogram oc- 
curs at the period 3.9078±0.0015 days (similarly observed 
in the NRS test), quite consistent with the orbital period 
determined by Casares et al. (f2005|) (3.90603±0.00017 
days) from radial velocity measurements of the stellar 
companion. The error in this measurement was estimated 
from Monte-Carlo-simulated time series containing a si- 
nusoid above a random background. The peak is highly 
significant, with a post-trial probability of less than 10~ 15 
that it results from a statistical fluctuation. This chance 
probability was estimated via Monte-Carlo simulation of 
random fluxes and also random re-sampling of fluxes 
(Fig. [3]). In Fig. [1] (middle panel) we also show the ef- 
fect of subtracting the orbital period, which removes nu- 
merous satellite peaks that are beat periods of the orbital 
period with the various gaps present in the HESS dataset 
(1-day, 28-day moon cycle, 365.25-day annual), that is, 
rational fractions of beat periods added to the orbital pe- 
riod. Fig. Q] (bottom panel) also includes results obtained 
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on the neighbouring VHE 7-ray source HESS J1825— 137 
(Aharonian et al. 12005a) . which is in the same field of view 
(FoV) as LS 5039 and therefore observed simultaneously. 
The HESS J1825-137 periodogram does not show statisti- 
cally significant peaks, demonstrating that the significant 
peak is genuinely associated with LS 5039. 

The ephemeris of Casares et al. (12005|) . determined 
from Doppler-shifted optical lines (observed in 2002 and 
2003), shows the binary makeup (Fig. g]) of LS 5039 
as comprising a compact object of mass > 1.38 M Q Q, 
in an eccentric e = 0.35 orbit around a stellar com- 
panion of mass ~20 Mq (with bolometric luminosity 
L* ~ 10 39 erg s^ 1 ). The separation (centre-to-centre) 
between these two components varies between 2.2i?» 
at periastron (<f> = 0.0 with reference epoch Tb(HJD- 
2400000.5)^51942.59) to 4.5i?* at apastron (0 = 0.5), for 
a stellar radius i?* = 7 x 10 11 cm. A range on the system 
inclination angle of 13° < i < 64° is inferred from the bi- 
nary mass function, the companion rotation velocity, the 
lack of X-ray eclipses (which assumes that the X-ray emis- 
sion occurs very close to the compact object) and lack of 
Roche lobe overflow. 

The phasogram (Fig. [5] top) of integral fluxes at en- 
ergies E > 1 TeV vs. orbital phase (</>) obtained on a 
run-by-run basis (one data run is ~28 minutes) clearly 
indicates that the bulk of the VHE 7-ray emission is con- 
fined to roughly half of the orbital period, covering the 
phase interval cf> ~0.45 to 0.9. The VHE flux maximum 
appears to lag somewhat behind the apastron epoch, and 
aligns better with inferior conjunction (0 = 0.716) of the 
compact object. Inferior conjunction of the compact ob- 
ject occurs when it is lined up along our line-of-sight in 
front of the stellar companion. The VHE flux minimum 
occurs at a phase (j> ~ 0.2, slightly further along the orbit 
than superior conjunction {(j) — 0.058), which is when the 
compact object is lined up behind the stellar companion. 
Note that the inclination upper limit i < 64° implies that 
direct views of both compact object and stellar compan- 
ion are always available. We define here two broad phase 
intervals for further study: INFC (0.45 < < 0.9) en- 
compassing inferior conjunction, and SUPC (<j> < 0.45 
and <j) > 0.9) likewise for superior conjunction. The phase 
error (A(f> = (TAP)/P 2 = 0.01) due to uncertainties 
in the period measurement AP = 0.00017 from Casares 
et al. PPDSI) . the dataset length T ~ 1000 days and 
P = 3.9 days, appears to be negligible. Nevertheless, fur- 
ther near-future optical line observations bracketing ours 
at VHE 7-ray observations would be desirable to check for 
the presence of systematic drifts in the orbital period. 

3.2. Phase-Resolved Energy Spectra 

The energy spectrum of the VHE 7-ray emission, and in 
particular how it might vary with orbital phase, is an 

1 A black hole of mass 3.7l}'g Mq for the compact object 
was derived by Casares et al. (|2005[> under the assumption of 
pseudo-synchronisation of the binary components 



important diagnostic tool. The differential photon energy 
spectrum (see Fig.il) (°- 2 to 10 -° TeV ) for INFC is con- 
sistent with a hard power-law where T = 1.85 ± 0.06 sta t ± 
0.1 S y S t with exponential cutoff at E = 8.7 ± 2.0 TeV (for 
fitted function dN/dE ~ E~ r exp(—E/E )). In contrast, 
the spectrum for SUPC is consistent with a relatively 
steep (r = 2.53 ± 0.07 s tat ± 0.1 sys t) pure power-law (0.2 
to 10 TeV) (see Fig. [6J. The spectra from these phase 
intervals are mutually incompatible, with the probability 
that the same spectral shape would fit both simultane- 
ously being ~2x 10~ 6 . Fitting a pure power-law (which 
is statistically sufficient at present) to narrower phase in- 
tervals of width Acj) = 0.1, and restricting the fit to en- 
ergies E < 5 TeV to reduce the effect of any cutoff, also 
demonstrates that a harder spectrum occurs when the flux 
is higher (Fig. [5] middle and bottom panels). Notably, the 
VHE flux at E ~ 0.2 TeV appears to be quite stable over 
phases and the strongest modulation occurs at a few TeV. 

We found no evidence for long-term secular variations 
in VHE flux on a yearly scale independent of the orbital 
modulation (Fig. [7]). The orbital modulation represents a 
VHE 7-ray luminosity (0.2 to 10 TeV; at 2.5 kpc) variation 
between 4 to 10 xlO 33 erg s _1 . Spectral fits and other nu- 
merical results are summarised in Tab.Q] Finally, in Fig. [8] 
we show that the spectral results obtained from a pure 
Hillas-based analysis using a separate calibration proce- 
dure are consistent with results from the Semi-analytical 
Model+Hillas results within statistical errors. 

4. Discussion 

The basic paradigm of VHE 7-ray production requires the 
presence of particles accelerated to multi-TeV energies and 
a target comprising photons and/or matter of sufficient 
density. In microquasars, particle acceleration could take 
place directly inside and along the jet, out to parsec-scale 
distances, and also at jet termination regions due to in- 
teraction with ambient matter (Heinz & Sunyaev [2002). 
A non-jet scenario based on acceleration in shocks created 
by the interaction of a pulsar wind with the wind of the 
stellar companion has also been suggested (Maraschi and 
TrevesQJHU Tavani et al. fTMTl Dubus[2006b]). The nature 
of the parent particles responsible for VHE 7-ray emission 
is under theoretical discussion, with both accelerated elec- 
tron (Aharonian & Atovan [l998[ Paredes et al. l2006p and 
hadron (Distefano et al. 120021 Romero et al. 2003} scenar- 
ios proposed. Observationally, electrons (eg. Corbel et al. 
[20021 Angelini et al. 12003)) and hadrons (Margon[1984]) are 
both known to be present inside jets. 

The orbital modulation in LS 5039, with a peak flux 
around inferior conjunction, minimum flux around supe- 
rior conjunction, and hardening of energy spectra, pro- 
vides new information about the physical processes in mi- 
croquasars. Our results can be compared with those 
at X-ray energies (3 30 keV), where interestingly, 
a spectral hardening with flux is also observed, as 
well as an indication for higher fluxes at <f> ~ 0.8 
in the phase-resolved light curve (Bosch-Ramon 
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Table 1. Photon energy spectra and luminosity (0.2 to 10 TeV) of the VHE 7-ray emission of LS 5039 for different 
orbital phase intervals. The broad phase intervals INFC and SUPC encompass the inferior and superior conjunction 
epochs. The orbital phase cj> is calculated from the ephemeris of Casares et al. (|2005|) . The best-fit function is indicated 
for each phase interval. The errors quoted are statistical with systematic errors in the 10 to 15% range. 



Orbital Phase Interval 




N 


r 




E 


Luminosity^ 




xlO -12 


ph cm _2 s -1 TeV~ 






[TeV] 


[erg s" 1 ] 


Fit function: 


dN/dE = NE~ r 


exp(-E/E ) 




INFC (0.45 < 4> < 0.9) 




2.28 ±0.10 


1.85 ±0.06 


8.7±2.0 


l.lxlO 34 


SUPC {<(> < 0.45) and </> > 0.9 




0.91 ±0.07 


2.53 ±0.07 




4.2xl0 33 


All Phases (time averaged) 




1.85 ±0.06 


2.06 ±0.05 


13.0±4.1 


7.8xl0 33 



Power-Law fit for energies E <5 TeV (narrow phase intervals) 



0.0-0.1 


0.84 ± 0.17 


2.62 ± 0.23 


0.1-0.2 


0.46 ± 0.21 


3.08 ± 0.47 


0.2-0.3 


0.81 ± 0.13 


2.46 ± 0.20 


0.3-0.4 


0.78 ± 0.18 


2.77 ± 0.60 


0.4-0.5 


1.88 ± 0.27 


2.32 ± 0.17 


0.5-0.6 


2.64 ± 0.20 


2.04 ± 0.10 


0.6-0.7 


2.20 ± 0.23 


1.88 ± 0.13 


0.7-0.8 


2.44 ± 0.20 


1.90 ± 0.10 


0.8-0.9 


2.96 ± 0.17 


1.94 ± 0.08 


0.9-1.0 


1.21 ± 0.25 


1.84 ± 0.25 



f At 2.5 kpc distance. 



et al. 2005). The VHE 7-ray modulation is an unam- 
biguous sign that periodic changes in the VHE 7-ray ab- 
sorption and/or production processes are occurring, and 
we discuss briefly how these could arise along with issues 
concerning the location and size of the 7-ray production 
region. 



4.1. Gamma-Ray Absorption 

VHE 7-rays produced close enough to the stellar com- 
panion will unavoidably suffer severe absorption via pair 
production (c + c~) on its intense optical photon field. The 
cross-section for pair production is dependent on the an- 
gle 9 between the VHE 7-ray and optical photons, and 
has an energy threshold varying as 1/(1 — cos 9). The 
level of absorption therefore depends on alignment be- 
tween the VHE 7-ray production region, the star, and 
the observer, leading to orbital modulation of the VHE 
7-ray flux (Protheroe and Stanev fl"987l Moskalenko ll995| 
Bottcher and Dermer CHH Dubus HUDSaJ) . For the orbital 
geometry of LS 5039, in which apastron and periastron 
are near the inferior and superior conjunction epochs, the 
absorption effect provides flux maxima at the phase of in- 
ferior conjunction, when 7 photons are emitted towards 
the observer parallel to stellar photons (cos 6— 1), and 
flux minima at superior conjunction (Dubus [2006ap . This 
picture agrees well with the observed phasogram (Fig. [5]), 
suggesting that absorption plays an important role. An 
additional key expectation is that the strongest absorp- 
tion, and hence modulation, should occur in the energy 
range E ~ 0.2 to 2 TeV (Dubus l20U6"a|) . However, the ob- 
servations show that the flux at ~ 0.2 TeV appears quite 



stable, suggesting that additional processes must be con- 
sidered to explain the spectral modulation. 

4.2. Gamma-Ray Production 

VHE 7-ray emission can be produced by accelerated elec- 
trons through the inverse-Compton (IC) scattering of stel- 
lar photons of the companion star, and/or from acceler- 
ated hadrons through their interaction with surrounding 
photons and particles. The efficiency of VHE 7-ray pro- 
duction under this scenario will peak around periastron 
(<t> = 0.0), reflecting the minimal separation between parti- 
cle acceleration sites and targets, and higher target photon 
densities. The observed phasogram is in contrast with this 
behaviour. However, important influences on the energy 
spectrum can arise from variations of the maximum en- 
ergy of accelerated electrons and dominance in the radia- 
tive processes (IC and/or synchrotron emission) by which 
they lose energy (known as cooling). The high tempera- 
ture of the companion star means that IC 7-ray produc- 
tion proceeds primarily in the deep Klein-Nishina regime 
(where the IC cross-section is sharply reduced compared 
to the Thompson regime) . The maximum electron energy 
B max is derived by equating the competing acceleration 
and dominant radiative loss timescales. The hard VHE 7- 
ray spectrum (r < 2.5) implies IC is the dominant cooling 
mechanism. Above some high energy e (discussed shortly 
below), synchrotron losses however are expected to take 
over. For dominant IC cooling we have (Aharonian et al. 
2006a) E max cx (B/w) 3 ' 3 , where B is the magnetic field 
and w is the target photon energy density. Since w oc 1 /d 
(for binary separation d), and assuming B oc 1/d, E max 
will increase by a factor ~10 from periastron to apastron, 
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leading to a spectral hardening around the apastron phase 
(for LS 5039, also near inferior conjunction). For magnetic 
fields expected within the binary orbit, B ~ 1 G, elec- 
trons with energy above e » Q[(B /G)(d/ R*)}^ 1 TeV^ will 
cool via synchrotron (X-ray) radiation in preference to IC 
cooling. Much stronger synchrotron losses will then pro- 
duce a steepening of the VHE 7-ray spectrum (Moderski 
et al. I2005P and thus the d and B dependence of this 
changeover energy e could also introduce a spectral hard- 
ening at apastron. Spectral hardening could also come 
from the angular dependence of the IC scattering cross- 
section (Khangulyan & Aharonian 2005), which peaks at 
small scattering angles (ie. around inferior conjunction). 
It is also feasible to consider production via accelerated 
hadrons. For considerably larger magnetic fields, for ex- 
ample at the base of the jet where B ~ 10 5 G could be 
expected, protons can be accelerated to multi-TeV 
energies. In such high magnetic fields, the strong 
cooling of electrons will also render any IC compo- 
nents negligible (see for example Aharonian et al. 
2006a) . Accelerated protons interacting with stellar wind 
particles and possibly X-ray photons H associated with an 
accretion disk, provide a plausible hadronic origin via the 
decay of neutral pions. 

4.3. Pair Cascades 

For magnetic fields B < 10 G (Aharonian et al. [2006a) . 
an extra complication arises from pair cascades initiated 
by absorption of first generation VHE 7-rays (Aharonian 
et al. I2006al Bednarek I2006|) . Cascades significantly in- 
crease the transparency of the source at ~ TeV energies, 
reducing absorption effects. Interestingly we do observe a 
VHE signal (excess significance 6.1a; 79 7-rays) in a phase 
interval <f> = 0.0 ± 0.05 where VHE 7-rays are unlikely to 
escape due to absorption (Dubus I2006a|) . A cascade, or 
at least an additional, possibly unmodulated component 
from another process may be present. 

4.4. 7-/?ay Production Region: Location & Size 

If the phasogram profile is the result of absorption, the 
7-ray production region, or at least a large fraction of it, 
should be embedded within the stellar photosphere de- 
fined where the absorption optical depth t 77 is >1. Such 
an optical depth occurs within ~1 AU of the stellar com- 
panion. However this size constraint comes with a caveat 
since absorption is unlikely to be the sole process present 
as argued earlier. An unmodulated component, possibly 
explaining the (j> = 0.0 ± 0.05 signal could for example 
arise if parts of the source lay outside the stellar photo- 
sphere, or is always situated in front of the stellar compan- 

2 The electron energy at which synchrotron and IC cooling 
times are equal. 

3 The X-ray disk luminosity, ~ 10 34 erg s , is marginally 
sufficient under this process to meet by itself that of the VHE 
7-ray emission. 



ion along the line of sight. The size constraint is consider- 
ably smaller than the jets/outflows observed out to a dis- 
tance of ^1000 AU from the binary system (Paredes et al. 
2002), and similar in size to the smaller-scale jets/outflows 
observed out to -2 AU (Paredes et al. I2000j) . The size 
constraint is also significantly smaller than the ~0.3 pc 
(28 arcsec at 2.5 kpc) upper limit (la) on the source ra- 
dius based on the HESS angular resolution. 

5. Conclusions 

In conclusion, new observations by HESS have established 
orbital modulation of the VHE 7-ray flux and energy spec- 
trum from the XRB LS 5039. The flux vs. orbital phase 
profile provides the first indication for 7-ray absorption 
within an astrophysical source, suggesting that a large 
part of the VHE 7-ray production region lies inside the 
photosphere (within ~1 AU) of the massive stellar com- 
panion. However, not all of the observed effects can be 
explained by absorption alone. Modulation of the energy 
spectrum could arise from changes in the maximum ener- 
gies of electrons responsible for the radiation, changes in 
the dominant radiative mechanism, and/or scattering an- 
gle dependence of the inverse-Compton scattering effect. 
A VHE 7-ray signal near <f> = 0.0 may arise from pair- 
cascades, or an unmodulated component produced out- 
side the photosphere. These observations provide key in- 
formation about the astrophysics associated with particle 
acceleration processes and subsequent VHE 7-ray produc- 
tion in XRBs. In particular, we are now able to begin to 
explore in detail the complex relationship between 7-ray 
absorption and production processes within these binary 
systems. 
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Fig. 1. Top: Lomb-Scargle (LS) periodogram of the VHE 
runwisc flux for LS 5039 (chance probability to obtain the 
LS power vs. frequency). Inset: zoom around the high- 
est peak (pre-trial probability ^10~ 20 ), which corresponds 
to a period of 3.9078±0.0015 days, compatible with the 
ephemeris value of 3.90603±0.00017 days (vertical red 
line at 0.2560 days -1 on the inset) from Casares et al. 
(2005|). The post-trial chance probability of the orbital 
period peak is found to be less than 10 -15 (see Fig. [3]). 
Middle: LS periodogram of the same data after subtrac- 
tion of a pure sinusoidal component at the orbital period 
of 3.90603 days. The orbital frequency peak has been re- 
moved as expected, as well as significant satellite peaks 
(see text). Bottom: LS periodogram of the HESS source 
HESS J1825— 137 observed simultaneously in the same 
field of view. The middle and bottom panel results are 
consistent with that expected of white noise over the range 
of frequencies sampled. The dotted lines correspond to a 
10~ 2 pre-trial chance probability. 
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Fig. 2. Top: Normalised Rayleigh Statistic (deJager 
1994) periodogram calculated from run-wise HESS fluxes 
for LS 5039. The middle and bottom panels depict 
the NRS after subtraction of the orbital period and for 
HESS J1825-137, respectively (As for the Lomb-Scargle 
test in Fig. [[]). 
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Fig. 3. Density function of the chance probability of the 
Lomb-Scargle power determined by Monte-Carlo and re- 
sampling methods after subtraction of the orbital period 
sinusoid. The expected exponential density functions are 
also indicated (solid lines). In the Monte-Carlo method, 
~ 10 6 random time series were generated to produce a dis- 
tribution of the highest Lomb-Scargle power. The highest 
power obtained after these iterations reaches just above 
20, well below the power of 62 obtained from the unshuf- 
fled lightcurve at the orbital period. An extrapolation of 
the curves is therefore used to estimate the chance prob- 
ability for powers above 20. 
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Fig. 4. The orbital geometry (Casares et al. 2005) viewed 
from directly above LS 5039. Shown are: phases ((f)) of 
minimum (periastron) and maximum (apastron) separa- 
tion between the two components; epochs of superior and 
inferior conjunctions of the compact object representing 
phases of co-aligment along our line-of-sight of the com- 
pact object and stellar companion. The orbit is actually 
inclined at an angle in the range 13° < i < 64° with re- 
spect to the view above. VHE 7-rays (straight black lines 
with arrows) can be absorbed by optical photons of energy 
hv £ , when their scattering angle 9 exceeds zero. 
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Fig. 6. Very high energy 7-ray spectra of LS 5039 for two 
broad orbital phase intervals (defined in the text): INFC 
0.45 < (j) < 0.9 (red circles); SUPC < 0.45 and <j> > 0.9 
(blue triangles). The shaded regions represent the la con- 
fidence bands on the fitted functions (Tab. [I]). Both spec- 
tra are mutually incompatible with the probability that 
the same spectral shape would fit both simultaneously be- 
ing ~ 2 x 10~ 6 . A clear spectral hardening in the region 
0.3 to ~ 20 TeV is noticed for the INFC phase interval. 
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Fig. 5. Top: Integral 7-ray flux (F > 1 TeV) lightcurve 
(phasogram) of LS 5039 from HESS data (2004 to 2005) 
on a run-by-run basis folded with the orbital ephemeris of 
Casares et al. p005p . Each run is ~28 minutes. Two full 
phase {(f)) periods are shown for clarity. The blue solid ar- 
rows depict periastron and apastron. The thin red dashed 
lines represent the superior and inferior conjunctions of 
the compact object, and the thick red dashed line de- 
picts the Lomb-Scargle Sine coefficients for the period 
giving the highest Lomb-Scargle power. This coefficient 
is subtracted from the light curve in Fig [1] middle panel. 
Middle: Fitted pure power-law photon index (for ener- 
gies 0.2 to 5 TeV) vs. phase interval of width Acf> = 0.1. 
Because of low statistics in each bin, more complicated 
functions such as a power-law with exponential cutoff pro- 
vide a no better than a pure power-law. Bottom: Power- 
law normalisation (at 1 TeV) vs. phase interval of width 
A0 = 0.1. 



2004 



H 




2005 



Fig. 7. Top: Integral flux (E > 1 TeV) vs. time (MJD) for 
LS 5039 on a run- by-run basis. Bottom: After subtrac- 
tion of the orbital period of 3.9063 days (this is achieved 
by subtraction of the Lomb-Scargle coefficients for the se- 
lected period). The average flux (dashed lines) for the 
post-subtracted light-curve is consistent with a steady 
source with a chance probability of 4xl0~ 2 . In addition, 
the average flux levels determined exclusively for 2004 and 
2005 data differ at a 2a level (statistical errors) but are 
consistent when factoring in systematic errors of ~10% on 
a run- by-run basis. 
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Fig. 8. Energy spectra of LS 5039 separated into the same 
broad phase intervals as for Fig. [SI INFC 0.45 < </> < 0.9 
(red circles); SUPC <f> < 0.45 and <p > 0.9 (blue triangles), 
comparing results from the semi-analytical Model+Hillas 
(filled markers) and Hillas-only (open markers) analy- 
ses (Aharonian et al. 12006b")) . respectively. The Hillas-only 
analysis has made use of an independent calibration chain. 



